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I. INTRODUCTION 

A single electron trapped in a Penning trap 111 is a 
unique quantum system in that it allows the measure- 
ment of fundamental physical constants with striking ac- 
curacy. Recently, the electron cyclotron degree of free- 
dom has been cooled to its ground state, where the elec- 
tron may stay for hours, and quantum jumps between 
adjacent Fock states have been observed ||. It is there- 
fore evident that the determination of the genuine (pos- 
sibly entangled) quantum state of the trapped electron is 
an important issue, with implications in the very founda- 
tions of physics, and in particular of quantum mechanics. 
After the pioneering work of Vogel and Risken ||| , sev- 
eral methods have been proposed in order to reconstruct 
the quantum state of light and matter Q], which range 
from quantum tomography [Q through quantum state 
endoscopy to Wigner function determination from 
outcome probabilities |7j]. Also, different techniques ||] 
have been proposed which allow to deal with entangled 
states. 

In fact, entanglement Q has been defined as one of 
the most puzzling features of quantum mechanics, and 
it is at the heart of quantum information processing. 
Some fascinating examples of the possibilities offered by 
sharing quantum entanglement are quantum teleporta- 
tion Jl0|,p ^| , quantum dense coding |Jr3|, entanglement 
swapping Jl 3|] , quantum cryptography , and quantum 
computation |15| . A striking achievement in this field has 
been the recent entanglement of four trapped ions Jul . 

In the present work we propose to reconstruct the full 
entangled state (combined cyclotron and spin state) of an 
electron in a Penning trap by using a modified version of 
quantum state tomography. Previous proposals fin\ need 
the a priori knowledge of the spin state and therefore arc 
not able to deal with entangled states. Our method, on 
the contrary, has the ability of measuring the full (entan- 
gled) pure state of the two relevant degrees of freedom 
of the electron. In order to reach this scope, our method 
takes advantage of the magnetic bottle configuration to 
perform simultaneous measurements of the cyclotron ex- 
citation and of the z component of the spin as a function 
of the phase of an applied driving electromagnetic field. 



The complete structure of the cyclotron-spin quantum 
state is then obtained with the help of a tomographic 
reconstruction from the measured data. 

The present paper is organized as follows: In Sec. || 
we outline the basic mode l of an electron trapped in a 
Penning trap, while in Sec. Ill we describe the main idea 
of our reconstruction procedure. In Sec. IV and ^ we 
concentrate on the measurement of the spin and on the 
tomographic reconstruction of the cyclotron states, re- 
spectively. We present the results of our numerical sim- 
ulations in Sec. VI, and conclude briefly in Sec. VII. 



II. THE BASIC MODEL 

Let us consider the motion of an electron trapped by 
the combination of a homogeneous magnetic field Bo 
along the positive z axis and an electrostatic quadrupole 
potential in the xy plane, which is known as a Penning 
trap jlj . The corresponding Hamiltonian can be written 
as 
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where A = (B Q y/2, — B Q x/2, 0) is the vector potential, 
c is the speed of light, d characterizes the dimensions of 
the trap, V is the electrostatic potential applied to its 
electrodes, and q the electron charge. 

The spatial part of the electronic wave function con- 
sists of three degrees of freedom, but neglecting the slow 
magnetron motion (whose characteristic frequency lies 
in the kHz region), here we only consider the axial and 
cyclotron motions, which are two harmonic oscillators ra- 
diating in the MHz and GHz regions, respectively. The 
spin dynamics results from the interaction between the 
magnetic moment of the electron and the magnetic field, 
so that the total quantum Hamiltonian is 

Htot = huj z {d\d z + 1/2) + huj c {d\d c + 1/2) + ^^a z . 

(2) 

In the previous expression we have introduced the lower- 
ing operator for the cyclotron motion 
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where [3 — muj c /2h and oj c — qB /mc is the resonance 
frequency associated to the cyclotron oscillation. For the 
axial motion we have 
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where us 2 — qV /md 2 . In the last term of Eq. (||), a z is 
the Pauli spin matrix and tu s = {g/2)uj c . 

The obtained Hamiltonian (Q) is then made of three in- 
dependent terms. Even though the only physical observ- 
able experimentally detectable is the axial momentum 
p z , in the following both the cyclotron and spin states 
will be reconstructed. Considering the eigenstates of & z 
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we can write the most general pure state of the trapped 
electron in the form 



I*} =cikMI T)+c 2 |^)| I) , 



(6) 



\tpi) and \ip2) being two unknown cyclotron states. The 
complex coefficients c\ and c<i, satisfying the normaliza- 
tion condition |ci| 2 + |c2p = 1, are also to be determined. 

The electronic state ( |) possesses two very interesting 
features: first, if the cyclotron states l^i) and j^) are 
macroscopically distinguishable, is a typical example 
of Schrodinger-cat state |is|| . Second, the full state of the 
trapped electron is an entangled state between the spin 
and cyclotronic degrees of freedom (unless = \ip2))- 
Introducing the total density operator R = (\&| as- 
sociated to the pure state |\E'), we can express the cor- 
responding total density matrix R in the basis of the 
eigenstates (jf]) of a z in the form 



R = 



(7) 



whose elements are operators. Its diagonal elements rep- 
resent the possible cyclotron states, while the off-diagonal 
ones are the quantum coherences and contain informa- 
tion about the quantum interference effects due to the 
entanglement between the spin and cyclotron degrees of 
freedom. 

It is also possible to give a phase-space description of 
the complete quantum state of the trapped electron by 
introducing the Wigner-function matrix ]l9| whose ele- 
ments are given by 



Wij(a) = (S t3 (a - a)) = Tr[i%(a - 
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In the previous expression the operator- valued delta func- 
tion S(a — a) is the Fourier transform of the displacement 
operator D(a) — exp(ao^ — a*d). 



III. MEASUREMENT SCHEME AND 
RECONSTRUCTION PROCEDURE 

The basic idea of our reconstruction procedure is very 
simple: Adding a particular inhomogeneous magnetic 
field — known as the "magnetic bottle" field [Q — to that 
already present in the trap, it is possible to perform a 
simultaneous measurement of both the spin and the cy- 
clotronic excitation numbers. Repeated measurements 
of this type allow us to recover the probability ampli- 
tudes associated to the two possible spin states and the 
cyclotron probability distribution P(n c ) in the Fock ba- 
sis. The reconstruction of the cyclotron density matrices 
Pa = iV'iKV'il (* = 1,2) in the Fock basis is then possible 
by employing a technique similar to the Photon Num- 
ber Tomography (PNT) |l7|]2C| l which exploits a phase- 
sensitive reference field that displaces in the phase space 
the particular state one wants to reconstruct |2l[| . 

In close analogy with the procedure described in 
Refs. HfH, the coupling between the different degrees of 
freedom in Eq. (|^) is obtained modifying the vector po- 
tential with the addition of the magnetic bottle field 0] 
so that A takes the form 
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Such a vector potential gives rise to an interaction term 
in the total Hamiltonian, 
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where the coupling constant n = qb/mc is directly related 
to the strength b of the magnetic bottle field. 

Eq. (11) describes the the fact that the axial angular 
frequency is affected both by the number of cyclotron 
excitations h c — d\d c and by the eigenvalue of & z . In 
terms of the lowering operator a z for the axial degree 
of freedom, the Hamiltonian that describes the interac- 
tion among the axial, cyclotron, and spin motions can be 
written as 
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(8) where the operator frequency Cl z is given by 



where i,j — 1,2 and Sij(a — a) is an operator in the 
product Hilbert space H = H cy c <8> ?4pin defined as 



m 







i 5. " 




-I) 





(13) 



2 



The operator Cl z is the modified axial frequency which 
can be experimentally measured jlj after the application 
of the inhomogeneous magnetic bottle field. What is ac- 
tually measured is an electric current (which is propor- 
tional to p z ) that gives the axial frequency shift 0% @- 
One immediately sees that the spectrum of Cl z is discrete: 
Since the electron g factor is slightly (but measurably 0] ) 
different from 2, Cl z assumes a different value for every 
pair of eigenvalues of n c and <j z . 



IV. SPIN MEASUREMENTS 

If one can perform a large set of measurements of £7 2 
in such a way that before each measurement the state 
|\&) is always prepared in the same way, it is possible 
to recover the probabilities P(T) and P(l) associated to 
the two possible eigenvalues of <j z , namely |ci| and \cz\- 
Recalling Eq. (0), we have 



f(T) = Tr cyc [(t |*)(*| T>] = l c i| 2 



(14a) 



P(i) = Tr cyc [(| |*)(*| I)] = |c 2 | 2 = 1 - P(T) . (14b) 

However, this kind of measurement does not allow to 
retrieve the relative phase 9 between the complex coef- 
ficients ci and C2 in the superposition (@). We can then 
add a time-dependent magnetic field ho[t) oscillating in 
the xy plane perpendicular to the trap axis p], i.e. 



bo(t) = bo cos(wt)x + 6q sin(wt)y 



(15) 



The resulting interaction Hamiltonian in the interaction 
picture is 



obtained applying the driving field ([L5J) for a time t = 
Tr/2u>R. We can now repeat the spin measurements just 
as we have described above in the case of the unknown 
initial state Soon after h»o(i) is switched off, the 

magnetic bottle field is applied again and the spin mea- 
surement is performed. Repeating this procedure over 
and over again (with the same unknown initial state) 
for a large number of times, it is possible to recover the 
probabilities P(T) and P(|) associated to the two spin 
eigenvalues for the state |*(t)) of Eq. (|s|). Without loss 
of generality, we can assume Ci £ R, c 2 = \c?\e lB , and 
(V'ilV'2) = re l P , which yield 



P(T) = i[l + 2r|d||o2|sin(0 + /3)] 
P(I) = ^[l-2r| Cl ||c 2 |sm(0 + /3)] 



(19a) 
(19b) 



It is important to note that the probabilities P(T) and 
P(l) can be experimentally sampled and that the mod- 
ulus r and the phase (3 of the scalar product (ipxfyi) can 
be both derived from the reconstruction of the cyclotron 
density matrices p\\ and P22, as we shall explain in the 
next section. Thus we are able to find the relative phase 
9 by simply inverting one of the two Eqs. (^), e.g. 
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The resulting 7r ambiguity in the arcsin function in the 
right hand side of Eq. (||(]) can be eliminated by choosing 
a second interaction time F and repeating the procedure 
above. 
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In the above equations, H mt (t) is the interaction Hamil- 
tonian in a frame rotating at the driving frequency w, 
while = gqbo/Zrnc is the Rabi frequency. 

The evolution of the state (0) subjected to the Hamil- 
tonian (nw in the resonant case u> — u s , yields the state 
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V. TOMOGRAPHIC RECONSTRUCTION OF 
THE CYCLOTRON STATES 

Let us consider again the state of Eq. (^|) : every time 
fl z (and therefore a z ) is measured, the total wave func- 
tion is projected onto | t)l ri c) or | |)|n c ), where \n c ) is 
a cyclotron Fock state. We then propose a tomographic 
reconstruction technique in which the state to be mea- 
sured is combined with a reference field whose complex 
amplitude is externally varied (as it is usually done in 
optical homodyne tomography |j|,|22|-Q) in order to dis- 
place the unknown density operator in the phase space (a 
technique very close to the PNT scheme p9| , p0| , p5| ) . In 
particular, we shall sample the cyclotron density matrix 
in the Fock basis by varying only the phase tp of the ref- 
erence field, leaving unaltered its modulus \a\ p7| , ^l| , p6[ . 

Following Ref. |l7j], immediately before the measure- 
ment of Cl z , we apply to the trap electrodes a driving 
field generated by the vector potential 
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where e is the field amplitude, which gives rise to a Hamil- 
tonian term of the form 



#drivc = -ift(ce- <u «*at 



(22) 



The time evolution of the projected density operator pa 
(i = 1,2) according to the Hamiltonian ( |22| ) may then be 
written in the cyclotron interaction picture as 



pu(e,t) = exp y-Hdiivctj p li (0)exp 



n 



(23) 



where we have defined the complex parameter a = et (t 
being the interaction time) and Pdrive is given by 



-ffdrivc = -ih(eal - e*a c ) . 



(24) 



The right-hand side of Eq. (|23|) is then the desired dis- 
placed density operator, where the displacement param- 
eter a is a function of both the strength e of the driving 
field and the interaction time t. Thus we can interpret 
the quantity 
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as the probability of finding the cyclotron state \tpi) with 
an excitation number n c after the application of the driv- 
ing field of amplitude e for a time t. Fixing a particular 
value of a, and measuring n c , it is then possible to re- 
cover the probability distribution ( ^5|) performing many 
identical experiments. 

Expanding the density operator pa in the Fock basis, 
and defining N c as an appropriate estimate of the maxi- 
mum number of cyclotronic excitations (cut-off), we have 
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The projection of the displaced number state \n c , a) onto 
the Fock state |m) can be obtained (generalizing the re- 
sult derived in Ref. (p7|) as 



(m|n, a) — (m\D(a)\n 
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x exp {i(m — n)[(p — TT0(n — to)]} 
xir^M 2 ), " (27) 

where 9(x) is the Heaviside function, L£J is the associ- 
ated Laguerre polynomial, and p = max{m,n}, v = 
min{m, n}. Inserting Eq. (p7l) into Eq. (|26]), we get 
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where p = max{fc,n} and v = min{fc,n}. 

Let us now consider, for a given value of \a\, PW(n, a) 
as a function of ip and calculate the coefficients of the 
Fourier expansion 



P^(n,\a\) = -L ^ tbpP® > (29) 



for s = 0, 1, 2, . . .. Combining Eqs. ( B8h and (gg), we get 
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where we have introduced the matrices 
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with p — max{m + s, n} and v = min{m + s, n}. 

We may now note that if the distribution pW (n, a) 
is measured for n G [0, N] with JV > iV c , then Eq. @ 
represents for each value of s a system of iV + 1 linear 
equations between the N+l measured quantities and the 
N c + 1 — s unknown density matrix elements. Therefore, 
in order to obtain the latter, we only need to invert the 
system 

JV 

(m+s\p u \m) = J2M% n (\ a \) pis,i) (nM) > (32) 

where the matrices M are given by M — (G T G)~~ 1 G T . It 
is possible to see that these matrices satisfy the relation 
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for m, w! = 0, 1, . . . , N c — s, . . ., which means that from 
the exact probabilities satisfying Eq. ( |30| ) the correct den- 
sity matrix elements are obtained. Furthermore, combin- 
ing Eqs. (^9|) and (|32|), we find 
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(34) 



which may be regarded as the formula for the direct 
sampling of the cyclotron density matrix. In particu- 
lar, Eq. ([54]) clearly shows that the determination of the 
cyclotron state requires only the value of (p to be varied. 
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We now only need to reconstruct the off-diagonal parts 
of the total density matrix (^), i.e. pxi — P21 — IV'iKV'al- 
This can easily be done under the assumption that the 
initial unknown electron state is pure, as in Eq. (|^). 
Then, we have 



(pl2)n,r 



i(P2 
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(35) 



Writing = ££L a n |n) and I $2) = £ 
Fock basis, we can obtain the desired coefficients from 
the recursive relation 
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where without loss of generality we can set a\ — 
[{pii)!.!] 1 ^ 2 € i?. A similar relation yields the coeffi- 
cients b n of 1^2)- Finally, the scalar product in Eq. j35| ) 
may be determined as 



N 



(37) 
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VI. SIMULATIONS AND RESULTS 

In this section we show the results of numerical Monte- 
Carlo simulations of the method presented above, which 
allow us to state that this technique may be quite ac- 
curate also in the experimental implementation. To ac- 
count for actual experimental conditions, we have con- 
sidered the effects of a non-unit quantum efficiency rj in 
the counting of cyclotronic excitations. When 77 < 1, 
the actually measured distribution P(k,a) is related to 
the ideal distribution P(n, a) by the binomial convolu- 
tion Mi 



n—k 



where 
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(39) 



Eq. (53) is then modified as 

N c -s 

V^(n,\a\) = E 2£l(|aU)<™ + 4S«H , (40) 



m=0 



where V^ s ' 1 ^ (n, \a\) is again defined according to Eq. (29) 
but now with V^(n,a) in place of P^'{n,a). In addi 
tion, Gn.m(\ct\,ri) is defined as 



sit(|aU) = X>^)4:l(M) 



(41) 



The matrices <7njL(|ct|, 77) can then be inverted in the way 
described in the previous section to obtain the matrices 
a |, 77) that can be used to reconstruct the density 

matrix. 

In the above reconstruction procedure, we have to con- 
sider two possible sources of errors associated to any 
actual measurement process. First, we have noticed a 
strong correlation between the statistical error in the sim- 
ulated reconstruction and the absolute value \a\ of the 
coherent field amplitude applied to drive the cyclotron 
motion. When \a\ is small, the density matrix elements 
near to the main diagonal are accurately reconstructed. 
Progressively increasing \ a\, the reconstruction of the off- 
diagonal elements becomes more accurate, while the el- 
ements close to the diagonal present significant fluctu- 
ations: to compensate for this effect, it is necessary to 
increase the number of measurements as much as possi- 
ble. Another source of error stems from the truncation 
of the reconstructed density matrix at the value N c in 
the Fock space. Neglecting the terms with n > N c — s 
in Eq. ( |30| ) causes a systematic error. This error can be 
reduced by increasing N c , which however may also cause 
an increase of the statistical error. This suggests that for 
a given number of measurement events there is an opti- 
mal value of N c for which the systematic error is reduced 
below the statistical error. 

In the following we will present two examples of appli- 
cation of the above method. We shall show the simulated 
tomographic reconstructions of an entangled electronic 
state of the type 



|*)= Cl | 7 )|T)+ C2 ^|7e l «)|l) 



(42) 



in which c\, C2, 0, 7, and £ are real parameters, and 
1 7) is a coherent state of the cyclotron oscillator. When 
£ 7^ 0, the state (j4|) is entangled, and when £ = tt (with 
I7I ^> 1) it represents the most "genuine" example of a 
Schrodinger-cat state JjJ. We shall display the results 
of our simulations both in terms of the density matrix 
of Ec^ (§ and of the Wi gner-function matrix JL9||29[ of 
Eq. (|8j). The latter can be derived from the density ma- 
trix elements (pij) n ,m through the relation pc| 

Wijfay) = -Y / (p ij ) n , n (-ire- 2 (* 2+ ^L n [4(x 2 + 1 ?)] 

7T * — ' 
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+ - E (PtfW-irv^PC* + iy ^ m ' n 
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fc=0 



In Figs. and | we show the numerical reconstruction 
of an asymmetric superposition of the type of Eq. p2] ) 
with ci = 1/2, c 2 = v/3/2, 9 = tt, £ = tt, 77 = 0.9 and 
7=1. In Fig. |l]we plot the results concerning the density 
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matrix, while in Fig. g those concerning the Wigner func- 
tion. In each figure, the true distributions are depicted 
on the left next to the corresponding reconstructed ones. 
Both the density matrices and the Wigner functions are 
well reconstructed. 

In Figs. |3| and ^ analogous results are shown for a sym- 
metric superposition state (42) with c\ = c-i = y/2/2, 
9 = 0, £ = 7T, T) = 0.9 and 7 = 1.5. Again, the re- 
construction is faithful. We would like to emphasize the 
particular shape of p 12 and W12 in both the examples 
above. It is due to the quantum interference given by 
the entanglement between the two degrees of freedom: in 
fact, in absence of entanglement (£ = 0) p\2 would just 
be a replica of the diagonal parts p\\ and p2i- In both 
the cases considered (Figs, [l] and ||) the off-diagonal den- 
sity matrix elements are real, due to the particular choice 
of the cyclotron states. The imaginary parts of the re- 
constructed density matrices turn out to be smaller than 
1(T 3 . 

We have performed a large number of simulations with 
different states and several values of the parameters, 
which confirm that the present method is quite stable 
and accurate. In addition, and for all the cases consid- 
ered, the values of the parameters ci, C2, and 8 (which 
obviously do not enter the plots of Figs. are very 
well recovered, with a relative error of the order of 10 -5 . 



VII. CONCLUSIONS 

In this paper we have proposed a technique suitable 
to reconstruct the (entangled) state of the cyclotron and 
spin degrees of freedom of an electron trapped in a Pen- 
ning trap. It is based on the magnetic bottle configu- 
ration, which allows simultaneous measurements of the 
spin component along the z axis and of the cyclotron 
excitation number. The cyclotron state is reconstructed 
with the use of a tomographic-like method, in which the 
phase of a reference driving field is varied. The numer- 
ical results based on Monte-Carlo simulations indicate 
that even in the case of a non-unit quantum efficiency 
the reconstructed density matrices and Wigner functions 
are almost identical to the ideal distributions. An exper- 
imental implementation of the proposed method might 
yield new insight in the foundations of quantum mechan- 
ics and allow further progress in the field of quantum 
information [Blfl. 
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FIG. 1. Simulated tomographic reconstruction of the density matrix for the state of Eq. ( |42| ) with ci = 1/2, C2 = V3/2, 
8 — 7r, £ = 7r, and 7 = 1. The quantum efficiency is 77 = 0.9 and 10 6 data per phase have been simulated. In this simulation 
an amplitude \a\ = 0.7 of the applied reference field has been used, a) (pn)n,m, b) {p22)n,m, and c) (pi2)n,m. Here and in the 
following figures the ideal distributions are displayed on the left, whereas the reconstructed ones are displayed on the right. 
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